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1. Introduction 

Over fifty years after Anderson's seminal paper on localization [Ij , the common wisdom 
is that in one dimension any disorder drives the phenomenon of Anderson localization. 
Even when the strength of disorder is arbitrarily small, all the eigenstates of the 
Hamiltonian get exponentially localized. The same picture applies to a large extent 
to the two-dimensional well. 

In a number of low-dimensional disordered systems the localization length ^{E) 
is however known to diverge at special energies. This phenomenon is usually the 
manifestation of some symmetry. A well-known example is that of tight-binding models 
with off-diagonal disorder, both in one dimension [21 [3l HJ [5l |6l [7] and on higher- 
dimensional bipartite lattices [H HI [TOl [11]. In this context the localization length 
becomes infinite at the center of the band {E — )■ 0). This phenomenon is due to 
particle- hole symmetry {E ^ —E). In the one-dimensional situation, it was already 
investigated in the pioneering work by Dyson [2J on disordered harmonic chains, while 
the interest in two-dimensional models has been revived by a connection with the integer 
Quantum Hall Effect [121 [13]. 

It is natural to wonder what are the dynamical consequences of the divergence of 
the localization length at the band center in this class of models. The usual setting 
consists in looking at the spreading of the wave-packet of a particle initially launched 
at the origin. A variety of rigorous general results, following the pioneering work by 
Guarneri [H], relate the growth of the moments of the position of the particle to 
dimensions of the spectral measure. The latter results however have very little predictive 
power in the present situation. 

The goal of this work is to obtain quantitative, albeit non-rigorous dynamical 
results, which can be checked against numerical data. Our main finding will be that the 
quantum motion exhibits a peculiar kind of anomalous diffusion which can be referred 
to as bi-fractality [15]. The disorder- averaged density profile Pn{t) of the particle at 
time t becomes critical in the long-time regime, in the sense that it falls off as a 
power of distance, as p„ = lim(_j.oo Pn(^) ~ l/\n\'^°~^^. As a consequence, the moments 
Mj(^) = d'^l'')* = J2n\''^\'^ Pn(t) of the position of the particle diverge with time t 
whenever q exceeds qo, while they have finite asymptotic values for q < q^. For the 
chain with off-diagonal disorder, investigated in Section 2, we obtain go = 2. 

Another example of interest is the random-dimer model [T6l [T71 [TBI [I9] . This is 
a one-dimensional tight-binding model with diagonal binary disorder, so that the site 
energies Vn = ±^ assume the same values on dimers, i.e., pairs of consecutive sites. The 
localization length is known to diverge at the special energies E = ±V when ^ < 1, 
in units where the hopping amplitude is one. We demonstrate in Section 3 that the 
dynamics of the random-dimer model fits within exactly the same framework, with 
go = 1/2. These findings confirm earlier heuristic results [16], namely M2{t) ~ t^^"^ in 
the generic situation {V < 1), while M2{t) ~ t in the critical situation {V = 1). In the 
first case, our predictions are in agreement with recent rigorous results [T9] . 
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Section 4 is devoted to a summary and discussion of our findings. We also argue wliy 
the same framework is expected to describe the quantum dynamics on two-dimensional 
bipartite lattices with off-diagonal disorder. 



2. The chain with off-diagonal disorder 

In this section we investigate the time-dependent discrete Schrodinger (tight-binding) 
equation on the one-dimensional chain with off-diagonal disorder: 

iV'n = Vn1pn+1 + (2-1) 

The hopping amplitudes Vn between sites n and n + 1 are modelled as i.i.d. (independent 
and identically distributed) random variables. We setQ 

Vn = e'", 7^ = 0, il = crl (2.2) 
2.1. Spectral properties revisited 

The time-independent (spectral) problem, defined by the eigenvalue equation 

Etpn = Vnlpn+l + (2.3) 

has been the subject of much activity since the pioneering work of Dyson [2]. The 
tight-binding equation (12. 3p has the peculiar feature that the even and odd sub-lattices 
decouple at zero energy. Hence the general solution at E = can be written down 
explicitly: 

ij2p = (-l)Pe^^-V^o, ^2p+i = (-l)^e^i-^^''+Vi, (2.4) 

where Bn is the following random walk (discrete Brownian motion) associated with the 
disorder: 

n-l 

5„= J2i-^r^m. (2.5) 

m=0 

The wavefunction (12. 4p exhibits unusual scaling properties [5]. We have = a'^n, and 

so the zero-energy wavefunction typically grows sub-exponentially as In ~ an^^"^. In 
other words, the effective localization length grows with distance as 

^1/2 

an) , (2.6) 

cr 

SO that in the thermodynamic limit, the localization length is infinite at zero energy. 

This peculiarity affects the behavior of various quantities near the band center. The 
Lyapunov exponent 7(-E) at energy E and the integrated density of states H{E) up to 
energy E exhibit the following Dyson singularities [21 HI [B] as E 0: 



1{E) ^ ^ 



2 



In E 



r2 



2(ln|E 



(2.7) 



X Our definition of cr^ differs from tfiat used in [B] by a factor 4. 
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The first result implies that the localization length ^{E) = l/'y{E) exhibits a logarithmic 
divergence as energy E approaches the band center, whereas the second one shows that 
states pile up near the band center, as the density of states strongly diverges as — )■ as 

This expression is at the verge of being non integrable. 

It will prove useful to first consider the finite-size scaling consequences of the Dyson 
singularity. On a finite chain made of an even number N of sites, with periodic boundary 
conditions, the spectrum manifests an exact particle-hole symmetry [9l |20] : if ipn is an 
eigenvector for the energy eigenvalue E, then {—l)^ipn is an eigenvector for the opposite 
energy {~E). Furthermore the expressions (12 ■4p show that E = can be an eigenvalue 
only if is a multiple of 4. This eigenvalue is then twofold degenerate. The condition 
on the disorder for this to occur is B]\r = 0, i.e., V1V3 . . . Vat-i = V2V4 . . . V^. 

We denote by Ei the energy gap of a finite system, defined as the smallest positive 
energy eigenvalue. The form f l2.7p of the integrated density of states suggests to set 
Xi = —In El. Typical values of the latter quantity can be expected to be such that 
H{Ei) - 1/2 ^ crV(2Xf) ~ We thus obtain the finite-size scaling law 

Xi = -\nEi^aN^/^Yi, (2.9) 

where Yi is a fluctuating quantity of order unity. The eigenstates at energy ±Ei are 
expected to be the most extended ones. The Dyson singularity (12.71) implies that their 
localization length scales as ^1 ~ Xi/cr^ ~ N^^'^/a. This result is in agreement with 
the estimate (12. 6p of the effective distance-dependent localization length ^(n) at zero 
energy, taken at the maximal distance n ^ N. 

Figure [1] shows a plot of Xi/a against N^^'^. The energy gap Ei is obtained by a 
numerical diagonalization of the Hamiltonian matrix on finite chains of sites, with 
periodic boundary conditions. Only values of multiples of 4 have been used. [Data 
for = 2 (mod 4) exhibit stronger corrections to the asymptotic linear behavior. This 
phenomenon can be related to the preceding observation that E = can only be an 
exact eigenvalue if A^ is a multiple of 4.] In numerical simulations hereinafter, the £.„ are 
drawn uniformly in the interval —W/2 < En < W/2, with W"^ = 12cr^. We shall report 
results for three values of the disorder strength: a = 0.5, 1, and 2. Data are averaged 
over 10^ samples for each value of A^ and a. The observed linear behavior corroborates 
the scaling law (12. 9p . The slope of a common least-square fit to all datasets (blue line) 
yields iT ~ 0.98. 

Figure [2] shows a histogram plot of the distribution of the reduced variable 
y = Yi/Yi = Xi/Xi for the largest system considered, i.e., A^ = 100. This distribution 
is observed to fall off very fast, both at small and large values of the variable, and to 
have a rather weak dependence on a. 
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Figure 1. Plot of Xi/a against N^^^, for three values of the disorder strength a. The 
blue line is the result of a common least-square fit to all datasets. 
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Figure 2. Histogram plot of the distribution of y = Yi/Yi = Xi/Xi for N = 100. 



2. 2. Dynamical properties 

Let us now return to the quantum dynamics of the particle described by the time- 
dependent equation fl2.ip . We assume that the particle starts at the origin: 

^n(0) = 5„,0. (2.10) 

Our basic observables will be the moments of the position of the particle at time t: 



M,(t) = (|n|^), = ^HV„(t), (2.11) 

n 

where 



Pn{t) = \MtW (2.12) 

is the disorder-averaged density profile at time t, whereas q is an arbitrary positive index. 
Hereinafter, the brackets (■ ■ ■)t denote an average over the quantum state at time t 
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described by the wavefunction ipnit), whereas the bar ■ ■ ■ denotes an ensemble average 
over the disorder, i.e., over the distribution of the random hopping amphtudes {Vn}- 
On a finite sample made of sites, the moments 

{\n\')t = j:\^\'\Ut)\' (2.13) 

n 

can be expressed in terms of the N energy eigenvalues E°' and of the corresponding (real 
and normalized) eigenvectors ip'^. Expanding the initial condition as 

^n(0) = 5.,o = EC^o, (2.14) 

a 

we find the wavefunction at arbitrary time t > 0: 

Mt) = Y.^-''^^'^:K- (2.15) 

a 

Combining fl2.13p and fl2.15p . we arrive at 

(l^n* = E e-^(^"-^')Vo° E C <■ (2.16) 

The long-time behavior of the above exact expression can be analyzed in the 
following heuristic way, along the line of thought of [15]. First, neglecting interference 
terms between different quantum states, we obtain 

(i^r)*^E(C)'Ei^r(cf- (2-17) 

a n 

Second, we model the eigenstates by considering that the probability density {ipnY is 
roughly uniform over the range \n — n'^j^\ < where we have introduced the center-of- 
mass co-ordinate n'^^ and the localization length of eigenstate a, defined as 

<n. = E^(C)', = E(^ - <m)'(C)'- (2.18) 

n n 

This picture implies that the eigenstates which contribute to the sum (12.171) are those 
such that |^^cml — ^"'^ Putting everything together, we are left with the following 
asymptotic long-time estimate for the disorder-averaged moments: 



M, = ^hmM,(t)~-E(e)^- (2.19) 

Finally, for a large system, the sum over discrete states can be replaced by an integral 
over the spectrum: 

-1 ^{E)''p{E)dE. (2.20) 

In the customary situation of the Anderson model with diagonal disorder, the 
moments of the position are thus predicted to saturate to finite limiting values of order 
Mg ~ ^Q, where for definiteness is the localization length at the band center. More 
quantitative predictions can be worked out, including universal pre-factors in the weak- 
disorder regime [T5]. 

The present situation of the Anderson model with off-diagonal disorder is radically 
different. Indeed, the Dyson singularities (12. 7p which affect the localization length and 
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the density of states near the band center cause the integral in fl2.20p to diverges for 
sufficiently large q. More precisely, taking X = — In 3> 1 as independent variable, 
the expressions (ETD yield ^{E) ^ X/a^ and p{E) dE = dH{E) ^ a^dX/X^, and 
so (I2:20|) reads 

X\« dX 



M, = {\n\^)-a^J ^. (2.21) 

The integral is divergent for q > 2. 

Let us first consider the stationary regime of very long times, and focus onto the 
integrand in fl2.2ip . not paying attention to the integration bounds for the time being. 
It is natural to identify the integration variable X with the position \n\ of the particle, 
according to |?t,| = Xja^ . We can thus read off from (12.211) that the asymptotic disorder- 
averaged density profile reached by the particle in the long-time regime, 

p„ = lim p„(t), (2.22) 

falls off as the power law 

'^"-^f^- (2-23) 

Let us move to the late stages of the dynamics of the quantum particle. As the 
latter starts from the origin, its wavefunction spreads progressively over larger and 
larger distances. For a large but finite time t, the power-law density profile (I2.23P is 
therefore expected to be cutoff at some scale n*(t). The growth law of the crossover 
scale n*(t) can be estimated as follows. First, as a consequence of the uncertainty 
principle, the finiteness of the observation time t induces a finite energy resolution 
£^*(t) ~ l/t. The largest observable localization length at time t therefore scales as 
the localization length at energy E^it) near the band center. Hence it diverges as 

i^{t) (lnE*(t))/(T2 ~ (lnt)/(T2. Then, the cutoff scale n^{t) can be estimated by 

equating the above ^^it) to the effective localization length ^{n) of (12. 6p . We thus obtain 

t 2.24 

The crossover scale n^it^ provides the appropriate cutoff to the integral in (I2.2ip . 
The disorder- averaged moments M^it) of the position can now be evaluated by means of 
the estimates (I2.23P and (I2.24p . We are thus left with the prediction that the moments 
grow with time on a logarithmic scale for g > 2, whereas they saturate to finite values 
for q <2: 



M,{t) 



(Int) 



2{g-2) 



^2(,-l) 



{q > 2), 



finite (g < 2). 



We have in particular 



(\ntf , , (lnt)4 

M3(t) ^ aa^^, M4(t) ^ a^^^. (2.26) 
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The above analysis has been tested against numerical data generated by directly 
simulating the time- dependent tight-binding equation f l2.ip . Figure [3] shows a log- 
log plot of the disorder- averaged density Pn{t) against time t, for a = 0.5 and even 
distances n ranging from to 12. The convergence to asymptotic densities p„ is observed 
to be slower and slower for larger distances. For each distance n, a black symbol 
marks the crossover time t^{n), defined operationally as the time where Pn{t) = Pn/^ 
(e ~ 2.718281). The oscillations which can be seen in po{t), and to some extent in P2{t), 
are reminiscent of those which are present in the absence of disorder, where the time- 
dependent wavefunction reads ipnit) = (— i)"J„(2t), with J„ being the Bessel functions. 
We mention for further reference that the latter wavefunction exhibits ballistic peaks 
around the positions n = ±2t, whose width and height respectively scale as t^^^ 
and t~^^^. More precisely, setting \n\ = 2t + t^^^z, one finds [15] 

i>„(t) ^ri/3Ai(^), (2.27) 

with Ai{z) being the Airy function. 


-2 

-4 

S -6 

a 
Cl 

a -8 

-10 
-12 

-14 

-2 2 4 6 8 

Int 

Figure 3. Log- log plot of the disorder-averaged density Pn{t) against time t, for 
a = 0.5 and various even distances n. Symbols: crossover times t*(n). 

Figure H] shows a log-log plot of cr^ times the asymptotic densities pn, against 
distance n. The data corroborate the prediction (12. 23 p . shown by the blue line with 
slope (—3). The amplitude of the law may have some weak residual dependence 
on a, besides the main scaling in 1/cr^. The corrections to the asymptotic behavior are 
observed to be larger at weaker disorder. 

Figure [5] shows a plot of {\nt^(n)) / a against n^/^, where the crossover time t^(n) 
has been defined by the condition p„(t) = Pn/e, as said above. This plot is a dynamical 
analogue of Figure [H where the system size is replaced by the distance n traveled by 
the particle, whereas time t replaces the inverse of the energy gap Ei. It is clearly far 
more difficult to get accurate dynamical data than spectral ones of the same quality. In 
spite of this, and although the available data span a rather limited range, they however 
point toward an asymptotic linear growth, whose slope is nearly independent of disorder. 
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Figure 4. Log-log plot of cr^ times the asymptotic densities pn, against distance n. 
The blue line with slope (—3) illustrates the theoretical prediction (|2.23p . 

These observations support the prediction (12.241) . again up to a possible weak residual 
dependence on a. The slope of the blue line has been chosen for definiteness to be the 
same as that of Figure [TJ 
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1/2 

n 

Figure 5. Plot of (lni,(n))/(T against n^^^. The slope of the blue line is the same as 
that of Figure [TJ 

Figure E] shows plots of (a^M^it))^/^ (left) and (a^M^it))^/'^ (right), against Int. 
The data exhibit linear growth laws, whose slopes seem independent of a, at least to a 
good approximation, thus corroborating (12. 26 p . The slopes of common least-square fits 
yield 03 ~ 12 and a4 ^ 6. 

Figure [7] shows a plot of a'^M2(t) against (lnt)/cr. This quantity is predicted 
in (I2.25P to grow logarithmically, with an amplitude independent of the disorder 
strength a. The available data cannot lead to a firm conclusion in this marginal 
situation. They however show a trend toward the predicted behavior, shown by the 
blue line (up to an unknown multiplicative factor). 
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Figure 6. Plot of (a^Msit))^^'^ (left) and {a^Mi{t)y/^ (right), against Int. Straight 
lines: common least-square fits yielding the amplitudes given in the text. 
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(In t)/a 

Figure 7. Plot of a'^M2{t) against {lnt)/a. The blue line illustrates the expected 
logarithmic growth (|2.25|) . 

The estimate fl2.24p of the dynamical crossover scale n^{t) provides the opportunity 
of putting the present problem in perspective with the Sinai walk [21j, namely the 
diffusion of a classical particle in a one-dimensional quenched random potential W{x), 
modeled as a Brownian motion such that W{xy = a'^x. The Sinai particle is known to 
have a very slow anomalous motion, such that x*(t) = i{x{ty))^^'^ ~ (lnt)^/cr^. This 
estimate bears a striking resemblance with our crossover scale n^,{t). Both models are 
indeed known to share a certain number of common features j22j . 

The different behaviors of the moments Mq{t) for q < qo and q > qo emphasized 
in fl2.25p (with go = 2) can be referred to as bi-fractality. This word was already used 
in the framework of a free quantum particle described by a tight-binding model in one 
dimension [15]. In this situation, the moments of the position of the particle exhibit 
ballistic scaling, i.e., Mg ~ for all positive q. The moments of the probability density 
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however exhibit a non-trivial scahng: 



^ (g = 2), (2.28) 

^-(2,-l)/3 > 2). 

This behavior is due to the balhstic peaks of the wavefunction tpnit) around n = ±2t 
[see (I2.27P ]. These peaks bring neghgible (resp. dominating) contributions for q < 2 
(resp. q > 2). The word bi- fractal had already appeared in the context of the distribution 
of matter in the Universe [23] . 

Coming back to our findings, the main consequence of the divergence of the 
localization length at the band center turns out to be the critical nature of the asymptotic 
density profile of the particle in the long-time regime, which falls off as a power of 
distance, as p„ ~ As long as we probe a moment Mq(t) with q < qo, the 

fat power-law tail of the latter profile is irrelevant, and so Mq{t) remains finite. The 
marginal moment (g = go) slowly diverges with time, as the logarithm of the cutoff 
scale n^{t). For q > qo, the divergence of the localization length is somewhat faster, as 

Finally, there is a qualitative difference between the kind of anomalous diffusion 
implied by the bi-fractal behavior met in the present work [see fl2.25p . fl3.16p . f l3.19p ] and 
conventional anomalous diffusion. The latter phenomenon is usually characterized by a 
single growing length scale i{t), so that Mg(t) ~ i{t)'^ for all positive q. In the case of 
self-similar growth, like e.g. in fractional diffusion, one has i{t) ~ t", and so Mq{t) ~ t°"^ 
for all positive q [see [21] and the references therein]. 

3. The random-dimer model 

We now turn to another example of a one- dimensional tight-binding model with unusual 
localization properties, namely the random-dimer model [161 [13 IHl [E]. The time- 
dependent Schrodinger equation reads 

^Ipn = + 1pn-l + K^n, (3.1) 

where the site energies Vn follow a sequence of dimers, i.e., they take random values on 
pairs of consecutive sites. In other words, 

V2k = V2k+i = Veu, (3.2) 
where takes values +1 and —1 with equal probabilities, and ^ > for definiteness. 

3.1. Spectral properties revisited 

The time-independent problem is defined by the eigenvalue equation 

Elpn = + + KV'n- (3.3) 

It was shown first by Dunlap et al ^16j that the random-dimer model exhibits a 
transparency phenomenon at E = ±V, and extended states near these special energies 
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for V < 1. This phenomenon is best described in terms of transfer matrices. The 
tight-binding equation (13. 3p can be recast as 



T(E,K)L:" , nE,V)= ^ \. (3.4) 






The generic locahzation properties of the customary Anderson model with diagonal 
disorder, where the site energies Vn are i.i.d. variables, can be related to the fact that 
the transfer matrices associated with single sites do not commute: 

[T{E,V,),T{E,V,)] = {V,-V,)I^^^ J). (3.5) 

The transfer matrices T{E, V)"^ and T{E, —Vy associated with the dimers however obey 

[TiE,Vf,TiE,-Vf] = 2ViE'-V')(^^^ . (3.6) 

This identity underlines that something special occurs for the dimer model ai E = ±V. 
Owing to symmetry it is sufficient to consider the case E = V. We have then 

T{E,Vf = -(^^^ (3.7) 

so that adding a dimer with e = +1 amounts to changing the sign of the wavefunction. 
The system therefore supports extended states at (and near) E = V, provided this 
energy lies in the spectrum of the other species of dimers {e = —1). For the sites 
belonging to the latter dimers, (13. 3 p reads 

{E + V)^„ = V-n+l + Vn-l- (3.8) 

For E = V, the dispersion law of the latter equation reads V = cos q. The condition for 
having extended states therefore reads ^ < 1. 

The following properties are known for the random-dimer model [16l [171 US [19] . 

• In the generic situation {V < 1), the special energy E = V is inside the 
band of (13.81) . The localization length is known to diverge quadratically and 
symmetrically as 

m - j^^, (E -> r), (3.9) 
whereas the density of states is regular at ii^ = V^. 

Neglecting localization effects in the immediate vicinity of the special energy, the 
dispersion law of (13.81) reads E + V = 2 cos q. As a consequence, the dynamics of 
the effectively extended states is characterized by the ballistic (group) velocity 

v = ^ =2(l-r2)i/2_ (3^^Q) 

E=V 

• In the critical situation {V = 1), the special energy E = 1 meets the band edge 
of (13. 8p . and so the velocity v vanishes. The properties of the spectrum are very 
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asymmetric in the vicinity of the special energy. The locahzation length and the 
density of states respectively diverge as 

m-Y^, P(^) - _^Ey/2 (E^l-), (3.11) 
whereas the density of states is exponentially small as -E — ?■ 1"*". 

Consider now, along the lines of Section 2, a finite sample made of an even number 
of sites, covered by N/2 dimers. The number of eigenstates which are effectively 
extended grows as A^^/^, both in the generic and in the critical situation [T7]. In the 
generic situation (V^ < 1), the eigenstates such that the difference 6 = E — V is smaller 
than A^~^/^ have a localization length ^ larger than N. As the density of states is finite at 
E = V, the number of such effectively extended states grows as NS ~ A^^/^. Different 
estimates yield the same outcome in the critical situation {V = 1). The effectively 
extended states are now such that the difference 6 = E — 1 is negative and of the order 
of 1/A^. As the density of states diverges according to ( 13. lip , the number of effectively 
extended states grows as 5"-'^/^ ~ A^^/^. 

The above property will be the cornerstone of our heuristic analysis of the dynamical 
problem. It is therefore worth checking its range of applicability by means of numerical 
data. Let E"' be the A^ energy eigenvalues of a finite chain made of an even number A^ 
of sites, and ip'^ be the corresponding (real and normalized) eigenvectors. We define the 
participation ratio of state number a as 

= ;F=r^ — 7. (3.12) 

n 

This quantity is known to provide a measure of the extension of the eigenstate a, i.e., 
of the number of sites which participate in the wavefunction ip", in the sense that ip'^ 
takes appreciable values [25]. The participation ratio scales as A^ for an extended state, 
whereas it is of order unity for a localized state. As a consequence, the disorder-averaged 
mean participation ratio 

provides an operational measure of the number of effectively extended states on a finite 
sample of size A^. 

Figure M shows a plot of P against N^^'^. Data have been obtained by means of 
a numerical diagonalization of the Hamiltonian matrix on finite chains with periodic 
boundary conditions. The observed asymptotic linear growth confirms that the number 
of extended states follows the predicted N^^"^ growth law, both in the generic and in the 
critical situation, with prefactors of order unity, down to system sizes as small as 10. 

3.2. Dynamical properties 



We now move to the investigation of the quantum dynamics of the random-dimer model. 
We again assume that the particle starts at the origin, according to fl2.10p . 
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Figure 8. Plot of the disorder-averaged mean participation ratio P against N^^^, for 
V — 0.6 (generic situation) and V — 1 (critical situation). Curves going almost exactly 
through the data points: fits to the form P = AiV^/^ + B + CTV^^/^. 



The generic situation (V < 1) 

The dynamical properties of the random-dimer model in the generic situation have been 
the subject of recent rigorous works [I9]. We proceed along an entirely different route 
and employ the following heuristic line of reasoning. 

On a large but finite sample of size A^, the model has a number of order N^^'^ 
of effectively extended states. The projection of the initial wavefunction fl2.10p onto 
extended states therefore scales as the weight of those states in the whole spectrum, 
i.e., as A^~^/^. This component of the wavefunction spreads ballistically with velocity 
w in a transient regime (t ^ A^). Now, considering an infinite system but a finite 
observation time t, we are led to the picture that some random projection of the initial 
wavefunction fl2.10p . carrying a fraction of order of the total intensity, spreads 

ballistically. This picture translates into the following scaling form for the disorder- 
averaged density profile: 

in the scaling regime where both n and t are large and proportional. 

The asymptotic disorder-averaged density profile p„ thus falls off as a power law: 

~ TZTy^ ■ (3-15) 



\n 



Figure [9] shows a log-log plot of the disorder- averaged density Pn(t) against distance n, 
for V = 0.6 (so that w = 1.6) and times t = 100, 200, 400, and 800. The data exhibit 
sharp ballistic peaks around the expected locations n = vt (arrowheads). Both the 
width and the structure of these peaks are reminiscent of the behavior f l2.27p of the 
wavefunction in the absence of disorder. For t > n/v, the data then soon saturate to 
the predicted asymptotic power-law profile (I3.15p . 




Figure 9. Log-log plot of the disorder-averaged density Pn{t) against distance n, for 
V = 0.6 and various observation times. Arrowheads: expected positions of the ballistic 
peaks (n = vt). Straight line with slope (-3/2): asymptotic power- law profile p.l5|) . 



The disorder-averaged moments Mq(t) of the position of the particle at time t grow 
with time for q > 1/2, whereas they saturate to finite vakies for q < 1/2: 

r (t;t)(2'?-l)/2 (g > 1/2), 

Mg{t) ~ I \n{vt) {q = 1/2), (3.16) 

[finite (g < 1/2). 

These predictions, and especially the power-law divergence with exponent (2g — l)/2, 
are in agreement with recent rigorous results [T9]. We thus have in particular 

Mi{t) ^ bit^/\ M^it) ^ h2t^'^. (3.17) 

Figure [lO] shows plots of Mi{t) (left) and M2{t)/t (right) against t^/^ for V = 0.6 
(so that V = 1.6) and V = 0.8 (so that v = 1.2). The plotted data exhibit linear 
growth laws with a very high accuracy. The corresponding amplitudes read bi ~ 4.2 
and 62 ~ 4.2 for V = 0.6, and 61 ^ 2.8 and 62 ~ 2.0 for V = 0.8. 



The critical situation (V = 1) 

Consider now the critical situation {V = 1), where the velocity v vanishes. It is worth 
mentioning that no rigorous result is available in this case. 

The previous line of reasoning can be adapted in the following heuristic way. 
On a large but finite sample of size A^, the model has a number of order A^^/^ of 
effectively extended states. The projection of the initial wavefunction (I2.10p onto those 
states therefore again scales as N^^^"^. The latter component of the wavefunction is 
again expected to spread ballistically in a transient regime, with a non-zero effective 
velocity v^{N) due to the finiteness of the system. The latter can be estimated as 
follows. The effectively extended states have an energy E such that 1 — E ~ 1/A^, 
whereas the velocity v is such that 1 — V ^ v'^/8 [see fl3.10p ]. Equating both estimates 
of the energy shift 1 — E 1 — V yields f (A^) ~ A^~^/^, implying therefore the existence 
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Figure 10. Plots of Mi{t) (left) and of M2{t)/t (right) against t'^^'^, for V = 0.6 and 
V = 0.8. Straight lines (slightly translated for readability): least-square fits yielding 
the amplitudes given in the text. 



of a crossover time t^,(A^) ~ N/v^:{N) ~ N^^"^. The growth law of this time scale is 
intermediate between the ballistic one {t ~ N) and the diffusive one (t ~ A^^). 

The above argument suggests the following scaling form for the disorder- averaged 
density profile: 

Pnit) ^ 7 ^1 = ^ ^1 • (3.18) 

The asymptotic disorder-averaged density profile p„ again falls off according to f l3.15p . 
The left panel of Figure [TT] shows a log- log plot of the disorder-averaged density Pn{t) 
against distance n, for times t = 100, 400, and 1600. The data converge toward the 
predicted asymptotic power-law profile fl3.15p . A comparison with Figure [9] however 
shows that the sharp ballistic peaks are replaced by broad shoulders. The right panel of 
Figure [TT] shows a rescaled plot of the same data. The very good quality of the observed 
data collapse, especially for the larger values of the ratio n/t^/^, strongly supports the 
scaling law (13.181) derived by heuristic means. 

As a consequence of (13.181) . the disorder- averaged moments Mq{t) of the position of 
the particle at time t grow with time for g > 1/2, whereas they saturate to finite values 
for g < 1/2: 

'^(2,-l)/3 (g>l/2), 

Mq{t) ~ J Int (g = 1/2), (3.19) 

[finite (g < 1/2). 

We have in particular 

Mi(t) ^ Clt^/^ M2(t) ^ C2t. (3.20) 

Figure fT2l shows plots of Mi(t) against t^/^ (left) and M2{t) against t (right). The plotted 
data clearly exhibit linear growth laws, thus corroborating (I3.20p . The corresponding 
amplitudes read ci ~ 3.2 and C2 ~ 5.3. 
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Figure 11. Left: Log-log plot of the disorder-averaged density Pn{t) against n, for 
V — I and various observation times. Straight line with slope (-3/2): asymptotic 
power- law profile p.l5|) . Right: Rescaled plot of the same data demonstrating the 
scaling law p.lSp . 
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Figure 12. Plots of Mi{t) against i^/^ (left) and of A'hit) against t (right) for 
V — 1. Straight lines (slightly translated for readability): least-square fits yielding 
the amplitudes given in the text. 



4. Discussion 

We have investigated two examples of disordered one-dimensional Hamiltonians which 
are exceptions to the common wisdom on Anderson localization, in the sense that the 
localization length diverges at some special energies. 

The main goal of the present study has been to provide a unified picture for the 
dynamical consequences of this divergence. The two examples we have considered 
are very different in nature. On the chain with off-diagonal disorder (Section 2), the 
localization length has a mild logarithmic divergence as — )■ 0. As a consequence, on 
a finite sample of N sites, the most extended eigenstate has an effective localization 
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length of the order of N^^'^/a. The system therefore supports no effectively extended 
state. On the contrary, in the random- dimer model (Section 3), the localization length 
has a strong power-law divergence as E ^ ±V. A large but finite sample now supports 
a large number of effectively extended states, growing as N^^'^. This difference probably 
explains why the study of the random-dimer model is far easier than that of the chain 
with off-diagonal disorder, both from the mathematical and the numerical standpoints. 

The main features of the quantum dynamics put forward in this work are 
nevertheless very similar in both models: 



The disorder-averaged densities pn = iV'nl^)!^ converge to non-trivial asymptotic 
values Pn, which fall off at large distances as a power law of the form 

[see (Km . (I37[5D . 

For a large but finite observation time, the power-law density profile is cutoff at 
some spatial scale n^,{t) which keeps growing with time. 

The disorder-averaged moments of the position of the particle keep growing with 
time as 

[\nn^{t) {q = qo), 

while they reach finite limits for q < qo [see (12.251) . (I3.16p . (13.191) ]. This kind of 
dichotomous behavior can be referred to as bi-fractality. 



Table 1. Summary of main results on the spectral and dynamical properties of a 
quantum particle on the chain with off-diagonal disorder and on the random-dimer 
model, both in the generic situation {V < 1) and in the critical one {V = 1). 



Model 


off-diagonal disorder 


generic dimers 
iV<l) 


critical dimers 
iV = l) 


Special energies Ec 





±v 


±1 


Localization length C,{E) 


-\n\E\/a^ 




1/\EtI\ 


Asymptotic profile p„ 




1/^3/2 


1/^3/2 


Index go 


2 


1/2 


1/2 


Crossover scale ?T,*(t) 


(lnt)Va2 


vt 


^2/3 


Moment Mg(t) (g > go) 


(lnt)2(5-2)/^2(,-i) 


(^t){2g-l)/2 


t(2g-i)/3 



Tabled] summarizes our quantitative findings concerning the chain with off-diagonal 
disorder and the random-dimer model. 

On higher-dimensional bipartite lattices, tight-binding models with off-diagonal 
disorder also exhibit specific features near their band center, again because of particle- 
hole symmetry [HI El [IHl HI]- In two dimensions, it has been argued that the energy 
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gap El on a finite sample of linear size L scales as Xi = —InEi ~ (InL)^. As a 
consequence, the integrated density of states exhibits an essential singularity of the form 

H{E) - 1/2 ~exp(-c(-ln|E|)^/^) (4.3) 

as E ^ 0, i.e., near the band center. The dynamical exponent has been first predicted 
by Gade [8] to be x = 2, while more recent studies [lOl [11] predict that the exact 
value of the dynamical exponent is somewhat smaller: x = 3/2. The dimensional 
scaling Ansatz H ~ 1/2 ~ 1/^^ suggests that the localization length diverges as 
^{E) ~ exp((c/2)(-ln|E|)i/^). These rough results fit within the same framework 
as our one-dimensional findings, with an asymptotic density profile falling off as 

a critical index go = 2, and bi-fractal moments scaling as 

rexp((g-2)(c/2)(lnt)2/3) (g > 2), 
M,(t)~ (lnt)2/3 (g = 2), (4.5) 

[finite (g < 2). 

One of the chief challenges for future work would be to establish by means of a more 
controlled approach the criticality of the stationary density profile and the validity 
of (14.51) in the two-dimensional situation. 
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Note added 



After completion of this article we became aware of the recent work [26] by Lepri et al. 
These authors investigate the spreading of a wavepacket in disordered harmonic and 
anharmonic chains. In the harmonic case they predict a power-law fall-off of the energy 
profile averaged over time and disorder, with a universal exponent equal to 5/2 for 
an initial displacement excitation, and to 3/2 for an initial momentum excitation. The 
anomalous spreading and the existence of a non-trivial limiting profile are demonstrated 
to be due to the divergence of the localization length of the phonon modes at low 
frequency. The work by Lepri et al and the present one therefore strongly corroborate 
each other, as they put forward essentially the same scenario. 

We thank Stefano Lepri for having made us aware of this very interesting piece of 
work. 
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